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Abstract 

We show that, if decaying gravitino dark matter is responsible for the PAMELA and ATIC/PPB- 
BETS anomalies in the cosmic-ray electron and positron fluxes, both a reheating temperature and 
a gluino mass are constrained from above. In particular, the gluino mass is likely within the reach 
of LHC, if the observed baryon asymmetry is explained by thermal leptogenesis scenario. 
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The PAMELA experiment [l[ reported an excess of the positron fraction above 10 GeV, 



which extends up to about 100 GeV. The excess could be a signa' 
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the annihilation or 
we consider 



me ixu 

IB, II, 



decay of dark mater. Among many decaymg dark matter mod* |3 
here the gravitino dark matter with broken .R-parity [2| (see also Refs. |6|). In fact, it was 
shown in Ref. that the gravitino decaying via the bilinear i?-parity violation can explain 
the PAMELA data. 

The positron spectrum needed to explain the PAMELA excess is rather hard. If the 
positron fraction continues to rise above 100 GeV, the cosmic-ray electron flux as well may be 
significantly modified at high energies. Interestingly enough, the ATIC balloon experiment 
collaboration js| has recently released the data, showing a clear excess in the total flux of 
electrons plus positrons peaked around 600 — 700 GeV, in consistent with the PPB-BETS 
observation [9|. It is highly suggestive of the same origin for the PAMELA and ATIC/PPB- 
BETS anomalies, if both are to be accounted for by dark matter. As will be shown in 
Appendix, the decaying gravitino scenario can actually account for both excesses. We focus 
on this scenario in this letter. 

The gravitino is assumed to be the lightest supersymmetric particle (LSP). In the presence 
of the i?-parity violation, its longevity is due to a combination of the Planck-suppressed 
interactions and a tiny .R-parity violating coupling. For the latter we assume the so-called 
bilinear i?-parity violating coupling, which is parametrized by a dimensionless coupling ^ 
defined as the ratio of the vacuum expectation values (VEVs) of the sneutrinos to that 
of the standard-model like Higgs boson, where the subindex i(= 1,2,3) denotes the flavor 
dependence, e, \i and r (see Ref. jlOj] for more details). We assume that the decay of an 
electron-type dominates over the others throughout this letter, i.e., K\ 3> ^2,^3, since one 
cannot fit well the sharp cut-off in the ATIC data otherwise. Then the mass and lifetime of 
the gravitino should be in the following range to account for the PAMELA and ATIC/PPB- 
BETS excesses: 

m 3/2 ~ (1.2-1.4)TeV, (1) 
r 3/2 ~ O (10 26 ) sec. (2) 

Since all the other supersymmetric (SUSY) particles must be heavier than the gravitino, 
we expect a typical mass scale for the SUSY particles, especially the gluino, may be out of 
the reach of LHC. This would be quite discouraging for those who expect SUSY discovery 



at LHC In this letter, however, we show that the gluino mass is bounded from above and 
is likely within the reach of LHC, if the baryon asymmetry is explained by the thermal 
leptogenesis scenario. 



Let us first discuss the gravitino production in the early universe 
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'he gravitino is 



produced by thermal scatterings,* 1 and the abundance is given by [13|, [14j, ll5j 

Y m ~ 4 x 10- £(T R ) In (-^-) + (r^) , (3) 



where #3 and M 3 are the SU(3)c gauge coupling and the gluino mass, respectively, and both 
are evaluated at a scale equal to the reheating temperature, T R , in Eq. ([3]). For simplicity we 
have dropped contributions involving the U(l)y and SU(2)l gauge interactions, which are 
subdominant unless the bino and wino masses, M\ and M2, are much larger than M3. Thus, 
the reheating temperature and the gluino mass are constrained from above for the gravitino 
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abundance not to exceed the observed dark matter abundance, VloMh 2 ~ 0.1143±0.0034 
In Fig.[T]we have shown the upper bound on the gluino mass and the reheating temperature, 
where we have included contributions from U(l)y and SU {2) L neglected in Eq. ([3]). We have 
imposed a requirement that the gravitino abundance should not exceed the 95% C.L. upper 



bound on the dark matter abundance. We used the code SuSpect2.41 |17j] to calculate 
the gravitino abundance and the physical spectra for the superparticles, with the following 
boundary conditions at the GUT scale ~ 2 x 10 16 GeV; tan/3 = 10, sgnj/i] > 0, the vanishing 
A-terms, the universal scalar mass m = 2 TeV for the squarks and sleptons, m 2 ^ = rn 2 Hd = 
5 x 10 5 GeV 2 , and the U(1) Y and SU(2) L gaugino masses M l = 3.5 TeV and M 2 = 1.8 TeV. 
Those parameters are chosen so that the gravitino is LSP * 2 

The origin of the baryon asymmetry is a big mystery of the modern cosmology. The 
thermal leptogenesis scenario [18] is appealing, and the reheating temperature must be 
higher than about 2 x 10 9 GeV Il9[ for the scenario to work. The precise value of the lower 
limit depends on flavor effects [20j and the mass spectrum of the right-handed neutrinos. 
The detailed study showed the lower bound as T R > 10 9 GeV, which is represented by the 
horizontal gray band in Fig. [D We can see from Fig. [T] that the gluino mass is bounded 



* The inflaton decay may also contribute to the gravitino abundance 

HQ 

We focus on the thermal 
production, since the non-thermal gravitino production depends on the inflation models. 
# 2 In the case of m 3 / 2 = 1.4 TeV, the gravitino is LSP for M3 > 600 GeV for the adopted parameters. This 
does not affect the following discussion. 



3 



10 



10 



> 



H 



io c 



to 

< 



Ml = 3.5TeV 
M2= 1.8TeV 



m 3 / 2 =1.2TeV 



1.4 x 10 9 GeV 



sj m 3/2 =1.4TeV 



500 1000 1500 2000 2500 

M3 [GeV] 



3000 



FIG. 1: The upper bounds on the gluino mass M3 and the reheating temperature Tr, for the 
gravitino mass m.3/2 = 1-2 TeV (solid red) and 1.4 TeV (dashed blue). The horizontal thick gray 
(thin orange) line shows the lower bound on Tr > 10 9 (1.4 x 10 9 ) GeV for the thermal leptogenesis 
to work. Here we set the bino and wino masses to be M\ = 3.5 TeV and M2 = 1.8 TeV. We 
also show the gluino mass in the low energy, m g i u i no = 1.5, 2.0, 2.5, 3.0 TeV, as the vertical dotted 
(green) lines from left to right . 

from above, M 3 < 1.5 TeV at the GUT scale, for Tr to satisfy the lower bound * 3 . This 
constraint can be translated into that the gluino mass should be lighter than about 3 TeV in 
the low energy, taking account of the renormalization group evolution. If we take a slightly 
tighter bound on T R , say, T R > 1.4 x 10 9 GeV, for which the leptogenesis becomes easier, the 
gluino mass in the low energy must be lighter than about 2 TeV for 1713/2 = 1.2 TeV. This 
is a surprising result. If the ATIC anomaly is to be explained by the decay of the gravitino 
dark matter, we may worry that the SUSY particles are so heavy that they may not be 
produced at LHC. However, if we believe in the thermal leptogenesis scenario and impose 
the lower bound T R > 1.4(1) x 10 9 GeV, the gluino mass turned out to be lighter than about 
2(3) TeV. This is a good new for those who anticipate the LHC to discover SUSY. 



# 3 The upper bound on the gluino mass was also discussed in Refs. 2l|, |22( in different contexts 
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Several comments are in order. In the presence of the i?-parity violation, it is quite non- 
trivial whether the SUSY particles can be detected at LHC, even if they are produced. If 
the gluino is the standard-model LSP, they will escape the detector before it decays. The 
collider signature will look like a split SUSY model 23| . and it is not easy to collect and 
analyze those collider data properly. On the other hand, if the lightest SUSY particle is the 
neutralino, we will observe a large missing transverse energy. Note that we cannot impose 
the GUT relation on the gaugino masses, Mi = Mi = Ms, since the bino would be lighter 
than the gravitino in the low energy. We have implicitly assumed that Mi and M 2 are 
not much larger than M3, throughout this letter. Our argument will not be significantly 
modified unless Mi and M 2 are much larger than M 3 . 

In order to realize the lifetime ([2]) the Ki must be chosen to be K\ ~ 10~ 10 . Such a tiny 
i?-parity violation can be realized in a scenario that the -R-parity violation is tied to the 
B — L breaking [241]. 

In this letter we have argued that, if the both PAMELA and ATIC/PPB-BETS anomalies 
are accounted for by the decaying gravitino dark matter, the gluino mass as well as the 
reheating temperature are bounded from above. In particular, the gluino is likely well 
within the reach of LHC if we assume the thermal leptogenesis scenario. Unexpected good 
news from the indirect dark matter experiments may be indicative of a bright future in the 
new physics search at LHC. 



APPENDIX A: THE DECAYING GRAVITINO AND THE PAMELA AND 
ATIC/PPB-BETS EXCESSES 

Here let us show that the decaying gravitino of mass 1.2 — 1.4 TeV can account for both 
the PAMELA and ATIC/PPB-BETS excesses. For simplicity we assume the isothermal 
distribution for dark matter profile, although our results are not sensitive to the dark matter 
profile. The electron and positron obey the following diffusion equation, 

V ■ [K{E, r) V/ e ] + ^ [b(E, r)f e ] + Q(E, f) = 0, (Al) 

where f e is the electron number density per unit kinetic energy, K(E, r) a diffusion coeffi- 
cient, b(E, r) the rate of energy loss, and Q(E, r) a source term of the electrons. The electron 
and positron fluxes are related to the number density by $ = (c/47r)/, where c is the speed 
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FIG. 2: The predicted positron fraction for the decaying gravitino dark matter, together with the 
PAMELA data. 

of light. The analytic solution of Eq. flAll) was given in Ref. j^. In the following analysis 



we fix 1713/2 = 1.2 TeV and T3/2 = 1.0 x 10 26 sec. See Ref. [26| for the values of the diffusion 
constant and the energy loss rate, and the details of the diffusion model parameters. 

The bilinear .R-parity violating operators depend on the lepton flavor. We consider the 
gravitino decay of the electron-type, that is, K\ ^> k 2 , k%. In Fig. [2] we show the positron frac- 
tion for the three different diffusion models, Ml, MED and M2, together with the PAMELA 
data. The MED and Ml models give a better fit to the data. Similarly we show the pre- 
dicted electron plus positron flux together with the ATIC data in Fig. [3] We have adopted 
the background for the primary electrons and the secondary electrons and positrons given 



in Ref. 



27 



281 ]. with a normalization factor kt, g = 0.75 for the primary electron flux. As can 
be seen from Figs. [2] and El the gravitino dark matter can nicely fit both the PAMELA and 
ATIC data. 

For the i?-parity violating operators of the \i— and r— type, the PAMELA data may be 
explained, while they give a very poor fit to the ATIC data. 
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FIG. 3: The predicted electron plus positron flux for the decaying gravitino dark matter, together 
with the ATIC data. 
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